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> ■ Abstract 

oo 

Q^ ■ We study the effect of finite size of hadrons on the quark hadron phase transition 
!>; 

(<~^ . and in particular on the location of the critical end point of such a phase tran- 

. sition. The corrections to the hadronic equation of state are incorporated in a 
thermodynamic consistent manner for a vander Walls like interaction. For quark 



^ ■ gluon plasma phase we take a bag model like equation of state which takes into 

account the perturbative interactions among the plasma constituents. We find 
that for finite sized baryons the first order quark hadron phase transition is not 
possible for the entire QCD phase diagram. The end point of first order phase 
coexistence line arises towards the higher chemical potential values in comparison 
to the point sized baryons, beyond which the transition from hadronic phase to 
the quark gluon plasma phase might be either crossover or second order phase 
transition. Our findings are consistent with the finite size scaling ( FSS ) analysis 
of RHIC data which negates the critical end point with baryon chemical potential 
values of less than 400 MeV 
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Introduction 



With the discovery of a possible quark-hadron phase transition under the extreme 
conditions of temperature and pressure P], a numerous efforts have been put forward 
to map the QCD phase diagram and in particular to locate the critical end point of 
such a phase transition, where one expects the nature of this transition to change from 
first to second order or to a smooth crossover [2]. However various studies done so far 
predict different results regarding the location of critical end point in such a deconfining 
phase transition. This variation stems directly from the treatment of the underlying 
physics close to the phase transition region involved in these studies, which on a more 
fundamental level requires a theoretical description in terms of not yet fully understood 
quantum chromodynamics (QCD) in its non-perturbative regime. 
Keeping this into account an equation of state for a strongly interacting hadronic matter 
was developed in [3] , where it was argued that for a strongly interacting hadronic matter 
nearing the point of phase transition, one can in principle ignore the contribution of one 
pion and one kaon exchanges among various baryons in comparison to the scalar and 
vector meson exchanges, therefore one could model the strongly interacting hadronic 
phase with an effective Lagrangian density of the form L = Cisb + '^k.k. + '^mr- The EOS 
for such a strongly interacting hadronic matter was used to develop the quark-hadron 
phase transition curve, where it was found that a first order phase transition is possible 
only upto a certain point, commonly known as critical end point CEP, for such a phase 
transition, with coordinates coinciding to the CEP as found in one of the variants of 
lattice gauge theory, LR2 [1] 

However it was found that the mesons (a;,p, 0) responsible for repulsive interactions 
among different hadrons were found to be net baryon density hb dependent, which 
would imply that for a system with low baryon-chemical potential /is and high tem- 
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perature values, one can in principle generate large number of particle -antiparticle 
pairs without any significant repulsive interactions among them o This is exactly the 
region where the CEP of quark-hadron phase transition curve is located. Therefore 
any rectification corresponding to repulsive interactions are expected to effect the lo- 
cation of CEP of such a phase transition. In this work we focus on this problem and 
construct a relativistic mean-field ( RMF ) model, which incorporates the finite size 
volume correction in a thermodynamic consistent manner. 

In this section we present an equation of state ( EOS ), for a strongly interacting 
hadronic phase and a weakly interacting quark gluon plasma phase used in this work 
and their respective definition. 



I. Hadronic Phase 

To derive the EOS for the hadronic phase consisting of finite sized hadrons we proceed 
as follows. First of all we take into consideration a system of Boltzmann particles with 
two particle interactions. The quantum mechanical generalization of such an EOS as 
applied to the hadronic matter will be written in a form similar to the one derived using 
Relativistic-Mean Field theory with a correction term corresponding to the finite size 
of hadrons. 



II. Boltzmann Particles with two Particle Interaction 



Let us consider a system of Boltzmann-particles, with two particle interactions, in 



^The effect would be more pronounced for the case of zero baryon-chemical potential, where there 
are no repulsive interactions among hadrons for any temperature values. 



thermal and chemical equilibrium at a temperature 'T'. A thermodynamic description 
of this system can be obtained using a grand-canonical partition function, [5] 



\nZG = V 



n 



+ Y.B,{T)n^ 



here 'V is the volume of the system, 'n' is the particle number density and 'i?i(T)' are 
the viral coefficients. Using a single particle partition function '^o', 



Zq = g 



(2vr) 



T 



(2) 



Eq.(l) reduces to. 



Inzo = ln(n) + -^B,{T)n'-' 



1=2 



(3) 



Now if we define the expansions. 



Pd{n,T) = TY,Bi{T)n\ u,i{n,T) = tJ^ -^B,iT)n'-' (4) 

% -L 



i=2 



i=2 



the pressure 'P' and number density 'n' takes the form |6] 



P(T, /i) =^ In Zg = Tn{T, /i) + P,i{n, T) 



(5) 
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n{T,fi) =zoexp ( j = 9 J I f I (6) 

where 'g' and '/i' are the degeneracy and chemical potential respectively. Now one can 
easily verify that the functions, Uci{n,T) and Pci{n,T) are related by [6] 



JUci{n,T) ^ dP,i{n,T) 
dn dn 



Here it is quite clear that with the knowledge of the Pci{n, T) and Uci{n, T), one can in 
principle calculate the pressure 'P' using Eq.(5), however because of Eq.(7), only one of 
the two functions must be known a priori. For illustration consider a vander-walls EOS, 



nT a 

1 — VqU V? 



where n, is the particle number density and Vq is the volume of each particle. Equating 
with Eq.(5) one gets 



P^,{n,T)=nT-^^ + -^ (9) 
1 — f on 



It is clear that the contribution corresponding to the finite volume of particles is con- 
tained in first term, denoting it by P^^ and the second term by P^, which takes into 
account the attractive interactions among the particles, one can calculate the function 

Uci{n,T) as. 



UJn,T) 



therefore with the knowledge of P^q one can directly calculate the corresponding func- 
tion or vice versa. The generalization of above results to quantum statistics can be 
readily obtained and are as follows, [6] 



PiT, /i) 



gT f d^k 



(27r)3 



In 



1 + a exp 



(fc^ + m^)^ + Uq- 



(27r)3(A;2+^2)| 



-fik) + PQ 



Q 



(2^ 



fik) 



(12) 



and the distribution function is given by f{k) = jexp ^ (fc^+'»^)^+^Q M j ^ ^1 

. Here 

a= +1, -1 for fermions and bosons respectively. The functions Uq and Pq are the quan- 
tum mechanical analogs of the classical functions Uci{n,T) and Pci{n,T) respectively. 
It can be readily verified that for the limit of a — )■ 0, the pressure 'P' and number 
density 'n' ( Eq.(ll), Eq.(12)) approach their classical limit (Eq.(5), Eq.(6)) provided 
one has following set of relations 



limUg = Uci{n, T); limPg = Pci{n, T) 

a— >0 a— >-0 



(13) 



III. Quantum Statistics of Interacting Hadronic Matter 



For an interacting hadronic matter, comprising of baryons, pions and Kaons along with 
their anti-particles, in chemical and thermal equilibrium, we write the total pressure 
'P' as 



PQS = P = P^;s + Pb,l+PQ (14) 
where, the baryonic contribution to the pressure is 



a J (27r)3 ^ \ T 



and the bosonic contribution is. 



a , (2^^M' + "^"P f 



+ 



X In < 1 + a exp 



The above expression for the total pressure 'P', follows directly from Eq.(ll), where 
the contribution from all particles (hadrons) has been taken into account. Here, me, 
nib ^-re the masses of the free baryons and bosons carrying a chemical potential, /ib, 
//b respectively. For a system of interacting particles and anti-particles one expects 
these variables to change and take some new effective values. Such a variation in 
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the present model is controlled by the functions, Uq and Uq respectively. Here, for 
baryons, we sum up over the entire baryon octet {N, A, E, S, A) and for bosons, we 
include Kaons, k'^ ,k~ ,kP and pions tt"*", 7r°, 7r~ only. 
Now , for a system of interacting hadrons, the function Uq, takes the form. 



Uq =Uq {vb, Vb, T) 

^UQ{n",nB,nb,n^'',n^\T) 



(15a) 
(15b) 



where, vb, vi, are the effective chemical potential's of baryons and the bosons, n is the 
total number density of hadrons, ub and n& are the net-baryon density and net-boson 
density. Also n^'',n^'^ are the scalar density of bosons and baryons respectively. 
For a system of non-interacting particles, as soon as the interactions (repulsive/ at- 
tractive) are are turned on, it is clear that the chemical potential of each particle must 
change, and as the equilibrium state, starts to set in, the effective chemical potential for 
each particle should approach an equilibrium value. Thus for a system of interacting 
hadrons, in the equilibrium state, the chemical potential for baryons and bosons must 
take some equilibrium value, say, vb and Vb respectively. Therefore we define Uq as 
a function of effective chemical potentials, vb and Vb, rather than chemical potential 
corresponding to the free baryons and bosons. Further keeping into account the func- 
tional dependence of the classical function Ud, Eq.(4), one can safely assume that the 
quantum function Uq, with a limiting value of, lima_^o Uq — Ud, must also be a function 
of number density as well, in general. Now, for the hadronic system under considera- 
tion, which consists of baryons, pions and kaons, it is quite clear that the variable T, 
controls the total number density of hadrons in the system , and also the effective 
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chemical potential's, vb and Vb, control the net-baryon density and net-boson density. 
In addition, these variables can also be thought of to give rise to scalar baryon density 
n^'^ and an analogous quantity for bosons n^'' as well. Therefore it is quite feasible to 
define Uq as in Eq.(15b). This seems to us, the most general form of the field Uq. 
Now, we write the function Uq as, [fl 

Uq = Ui (n^, T) + U2 {n^^,n^\T) + U3 [ub, rib, n^\T) (16a) 

with the symmetry property, U3 — )■ —U3, under the transformation, Hb ^ —fJ'B or 

Therefore we define, 

Uq = Ui (n^, T) + U2 (n^^, r^^^ T) - U3 (n^, n^, r^^^ T) . (16b) 
Now for the thermodynamic consistency, it is clear that the thermodynamic relation 



dP 



B 



= riB (17) 



must always hold. Therefore for pressure 'P' as defined in Eq.(14), with the functions 
Uq and Uq as defined in Eq.(16a) and Eq.(16b), we haveQ 



^Appendix A.l 

*this follows directly from the symmetry of pressure under this transformation. 

5 here we have used, §^ = + + where, = Pg(n^, n^S n,, T) 



dP 
djiB 



Ub + 



dn 



H 



+ 



dPq 
"5 

OUb 



+ 



dP, 



Q 



— n 



H 



dn^ J d^B \dn^B dn^B j 



dUB 



dU^ \ dn 



- rib- 



duB J d/lB 



(18) 



Therefore for thermodynamic consistency one can write 



/OPq ^HdUi\ dn'' ^ / dPg dU2 \ dn^^ ^ / dPg dU^ ^ dU^ ^ Q 

\(?n^ dn'' J d^B \dn^B dn^^ J Q^^ ydns dnB dnB J diiB 



Now because of the independence of the variables {n" ,n^'^ ,nB) and therefore of the 
their derivatives with respect to /ib as well, the terms in the brackets should vanish 

separately. Equating first term to zero, we get, 

dn" dn" 

which after integration yields, 

Pq - Pqiin", T) + PqsK, n^, n^^, T) (19) 



Therefore the total pressure 'P' can be written as, 

Pqs = Pb;b + n,6 + Pqi + Pq2 (20) 



in 



where, the baryonic contribution P^;^ is, 




and for the bosonic contribution P^^^ we have. 




here E*^ — Eb + U2 and E^ — Ei, + U2 are taken to be the, in-medium effective energies 
of baryons and bosons. Also vb = IJ'B — U3 and Vb = fJ^b — Us represent the effective 
chemical potential of baryons and bosons respectively. 

IV. RMF Model and the Finite Size Correction 

In an RMF model for strongly interacting hadrons ( baryons + pions-l- kaons ),the 
Lagrangian density is 

C^otai ^jr^^jr^^jr^ (23) 



11 



where 



B 

- \S,,S^'^ + \ml<P,r (24) 



here S^''' = - R'^''' = d^'p'"' - d^p""^ + g pe"^" p^^" p""" and W^""" = df'u'' - & 

and the sum is carried over the entire baryon octet. 
Also 



^K=Y1 ^i^K^d^K - mlK^K - g^KrriKK^Ka - g^.KrriKK^Ka* - g^KKUXKu^ 

K 



gpKK^Tii dp Kpi^" - g^x^U dp Kcj)'' + {g^K^p + gpR^ipip + g^px^'i^fK'' K (25) 



and 



=229fi'^^d''n - mln^n - g^^m^n^no- - g^^nUdpTru'' - g pT.'n'^ r^i d pix pi^ 

TT 

+ {gui-K^^pi + g^-KTiPipf 7r"^7r (26) 
The equation of state EOS turns out to be [3] 



+ P{nB,rn,,n'^\n^^ ,T) (27) 
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The first two terms are the kinetic terms while the last term has been written in a 
generalized manner to show the dependency on baryon and boson densities. On com- 
paring with Eq.(20) it is quite clear that in the limit of vanishing hadronic volume, the 
two expressions become equal, i.e for vq — > 0, Prmft — Pqs, therefore one can readily 
correct the above EOS for finite size effects and the resulting expression becomes. 



Pi 



RMF 



+ P(r^B,n6,n^^r^^«,^) + P(r^^,^) (28) 
here the modified distribution functions for baryons and bosons /^^ are given by. 



/b(b) = S exp I '- I + o S (29a) 



■(fc2+mf)i+C/(n^,T)-M,(5) 



-1 



fh{b) = <^ exp — ^ \+a) (29b) 



with the effective chemical potential, vb — I^^b — QuB^ — 9(I>b4' ~ SpB^sbP, Vf, — /ii, — 
gojb^^ — g<i>b4> ~ gpbTsbP and the effective in- medium mass, m% = tub + gaBC + ga^cr*, 
ml — {ml + mb{gabO' + ga*b(^*))^ for baryons and bosons respectively. Here a, a*, uj, p 
and (f) are the mean-fields present in the system. The field equations can be deduced 
from the Lagrangian density Eq.(23) by minimizing the corresponding action and are 



as follows. For sigma (a) field 



m 



* 2 



m% {ub + ub) 



b={K,7r) 



for omega (u) field 



b={K,7T) 



+ CsUj'^J2g^B ^ J d^kiriB-riB) + Yl '^9cob 7^/"^^ 



X {E+nb + E-fib)] - Yl {'^^^9lb + '^9ojb9pbpr3 + 2g^b9m^) 

b={K,n) 



X 



lb 



{2t:Y J 2ub 



{rib + nb) 



for rho (p) field 



lb f d^k 



P=Y. 9pBTs / d'k {UB - nB)\ + 2g,,r, j 



X {E^Ub + E'rib)] - ^ i^g^^gpb^Tri + Ig^p^^rlp + 2g^b9pbT-i^) 

b={K,n) 

lb f d^k 



X 



{nb + Ub) 



for sigmastar (a*) field we have 



IB 

{2ny 



d^k 



m 



B 



m 



* 2 



(ub + Hb) 



- Yi 9a*bmb 

b={K,'K) 



X 



lb 



{2ny 



dH: , 

— {ub + nb) 
Zujb 
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and for phi (</>) field we have 



m, 



g<s>B 



B 
X 

X 



IB 

{2ny 



d^k {riB — tib) 



b={K,-K) 



lb 



_(2Tif J 2ujb 



{Eb'^nb + Eb rib)] - ^ {2g^bg'S>b(^ + ^g^bgpbT-sP + 2^glf^) 



b={K,Tr) 



lb 



(2vr) 



[ d^k . 

7 2^^"^ + 



(34) 



For the evaluation for pressure due to hadrons (pions + kaons + pions) all five coupled 
field equations are to be solved with two given parameters, temperature T and chemical 
potential /i. The coupling constant parameters are taken to be the one which have been 
deduced from low energy experiments. Here we want to make clear that it is not known 
what are the antiparticle couplings, we are assuming that it is of same magnitude as 
that of particles. 



V. Quark Gluon Plasma Phase 

For the weakly interacting QGP phase, we use a Bag-model [7] inspired EOS and restrict 
the description to three light quark-fiavors (u,d,s) and gluons with the perturbative 
corrections of the order of [8J. The pressure and energy density for this system are 
given by [9]. 



nj poo _ 

eQGP{T, =(1 - C) E ^ / ^'^^ V^fc^T^ [/.,/(A;) + f.jik)] 



4 A Ar„ 



oo 



+ (l-5C)^/ k'dkjk^ + mlf,{k)+ B (36) 



Here C = is the model parameter and accounts for the strength of interaction present 
in the system. Ng — 2{N^ — l),Nc are the number of transverse gluons and quark colors 
respectively, is the quark mas of flavor 'f ' and rUg is the gluon mass. The distribution 
functions for the quarks and gluons are fqj{k) and fg{k) respectively. Here 'B' is the 
bag constant and is second free parameter in our model. Now because of the energy 
scales present in the system , it is clear that one can treat two light quarks u and d 
to be almost mass-less. Therefore solving the above set of equations for two mass-less 
quarks (u, d) and an s-quark of finite mass m^, the expressions for pressure and energy 
density can be written as, 

PiT, ,) + N, + \,^T^ + ^,^) + i-^ £ dE iE^ - m.) ^ 

X (/fe + h)-B (37) 



£(T,/x) =3 



^ > +3iV 



90 



/ 



1-C 



dE 



xE'{E'-ml)'^ X {f, + f^)+B 



(38) 



Ifi 



Here Ng = 16 (l — |C) and Nj = 2 {1 — () are the effective number of gluons and the 
effective number of hght flavors. The third term gives the contribution of strange quark 
and strange anti-quark respectively. It is worth mentioning here that in this calculation 
we have taken two light quarks u and d to be mass-less and only strange quark to be 
of finite mass. 

VI. Results and Discussions 

To have a theoretical prediction about the QCD phase diagram we will now use the 
equation of state for the hadronic phase, corrected systematically to take into account 
the finite size of hadrons as explained above, together with the equation of state for 
the quark gluon plasma QGP. For hadronic system all five coupled field equations for 
the u,a,a*,p and field are to be solved in a self consistent manner, for two given 
parameters , temperature T and baryon chemical potential /i. For a given pair of (T, fi) 
values a complete set of field values is required to be obtained. Then using these field 
values different observable's of the hadronic system can be readily computed. In case 
of QGP phase the problem of obtaining the pressure and energy density is relatively 
simple and can be achieved quite readily. 

First of all we present some of the features of QGP phase used in this model. It is clear 
from the QGP EOS, that there are two model parameters, coupling constant ( = 
and Bag-value 'B', which must be determined uniquely to make any sensible prediction 
using this model. We fix these parameters by referring to some of the earlier work done 
in this direction. Following Ref. [31 [9] , we take the value of parameter zeta to be, 
( = 0.2, as it was shown that it is only for this value of zeta that the above defined 
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QGP equation of state shows a behavior consistent with the EOS derived using lattice 
QCD. Similarly for the parameter B we choose a value of B=344 MeV/ fm? following 
Ref. [9], where it was used successfully for the fluid-dynamical calculation of heavy- ion 
collisions. Here as already mentioned we will be treating light quark flavors (u , d) as 
mass-less and the quark 's' to be of finite mass. Using these parameters we first of all 
plot variation of Pqgp with interaction parameter C, as shown in Fig. 1. 




-400 I ^ ^ ^ ^ ^ 1 

50 100 150 200 250 300 

T(MeV) 

Figure 1: Pqgp v/s T, here C = 0.0, 0.2 for a and b respectively. Here /Xg=0 MeV. 

Here we show the variation of pressure in QGP phase for two different cases, a com- 
pletely non-interacting QGP phase {Q = 0.0) and a QGP phase with perturbative 
interactions {( = 0.2). As expected the pressure for the QGP system drops as the 
interactions are switched on. Starting from a negative value the pressure for the QGP 
phase becomes positive only after some finite temperature 'T', which is the general be- 
havior of bag model equation of state. It is in this region the vacuum pressure is larger 
than the pressure exerted due to quarks and gluons, which results in the instability of 
QGP. It is only when Pqgp > that the QGP phase is in a stable state. 
In Fig. 2 we show the variation of Pqgp with temperature 'T' for different values of 
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the quark-chemical potential fiq. Here we have kept the interaction parameter fixed at 
C = 0.2. It is this value of the parameter ( that we will use throughout our calculation. 
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Figure 2: Pqgp v/s T, C = 0.2 and fiq= 200, 170, 140, 120 MeV for a,b,c and d 
respectively. 



Next, we present some of the features of RMF equation of state. The values of various 
coupling constants used here are listed in the end. First of all we fix the function 
U{n^,T). From Eq.(lO) and Eq.(9) we can write, ( after neglecting the vander- walls 
type attractive interactions) 3 



U{n'',T)=T ^ """^^^ -Tln(l-t;on^) (39) 



1 — VqU 

where the corresponding function P{n^ , T) is 



P{n^T) = n"T-^^ (40) 
1 — Von" 



^In view of the strength of strong interactions one can easily neglect the contribution from such 
terms. 

1Q 



Here n is the total number density of hadrons and vq is the volume of each hadron. 
It is for simplicity that we have taken all hadrons to be of equal size. 
In Fig. 3 we show the variation of U{n^ ,T) with temperature for different values of 
baryon-chemical potential fiB- 




60 80 100 120 140 160 180 200 



T (MeV) 

Figure 3: Uydw v/s T, here fig = 190, 140 MeV for a and b respectively. 

Here we have taken only baryons to be of some finite volume Vq and for bosons we 
choose Vq = 0. This is consistent with pauli's exclusion principle from which it follows 
directly that it is possible for two or more bosons to occupy same position at the same 
time. 

In Fig. 4 we plot the variation of corresponding mean-field P{n^ , T) with temperature 
for different values of chemical potential. Among the possible values of baryon radius 
which lie in the range, = 0.5 — 0.8 fm [10], we have taken baryon radius to be 
0.8 fm and will use same value throughout our calculation. 

It is quite evident that the functions U^dw and Pvdw have a significant value at higher 

9n 




120 140 160 



T (MeV) 

Figure 4: P^dw v/s T, here Hq — 190, 140 MeV for a and b respectively. 

temperatures, therefore one should expect that the finite size correction should play an 
important role in the description of quark-hadron phase diagram. Next in Fig. 5 we 
show the variation of pressure in the hadronic phase with temperature for a fixed value 
of baryon-chemical potential /ib- 

In this figure we have compared the pressure for two different cases. In one we have 
considered baryons as point particles while in second we have taken baryons to be of 
finite volume vq. The drop in the pressure for a system of finite sized baryons can be 
attributed to the fact that for any non-negative finite value of mean- field U{n^ , T) the 

f I f I 

ratio of distribution functions, , ^ = e~^^ and . , ° = e~l^'-' are less than one, 

/b|i>()=o h\vQ=o ' 

P I 

which directly imply that, Pj^^JJ\ \ < 1- In the similar fashion one can also show 
that the P-T curves for other values of baryon-chemical potential as well. 
Now to develop the quark-hadron phase transition diagram we use the Gibbs-criteria of 
the phase-transition , which requires following set of relations to be valid at the phase 
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50 60 70 80 90 100 
T (Mev) 



Figure 5: P=Pb.mft v/s T, for fig =170 MeV. here hard core radius is tq = 0.8/m, 0/m 
for 'a' and 'b' respectively. 

equihbrium points, [TT] 

Pq = Ph, I^q = f^H, Tq = Th 

In Fig. 6 and Fig. 7 we show some of the phase co-existence points. The hadronic 
phase has been amputated beyond pt. B as the entropy density of hadronic matter 
for this region is negative and is hence unstable. The stability might arise due to the 
breaking of this entire hadronic phase into smaller components with pressure P > 0. As 
the temperature is increased further the system enters in the state with negative total 
pressure but with positive entropy density. This state is the hquid phase of hadronic 
matter which is appearing at very high temperatures. This observation for hadronic 
matter has been also reported in ref |12l |13]. For the baryon-chemical potential values 
less than 420 MeV the point of intersection starts to appear below the Prmft = 
axis. As for this region the pressure due to quarks and gluons is still less than the bag 



9.9. 



pressure 'B', the QGP phase does not exist. An increase of temperature by infinitesimal 
amount leads to the formation of QGP phase as can be seen from the large slope of 
P-T curves for this phase. The transition from hadronic to QGP phase therefore occurs 
without the appearance of the phase coexistence region and hence is a crossover. 
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120 130 140 
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150 160 



Figure 6: Prmft, Pqgp intersection plots. A, B, C are the intersection points. Here 
jiq— 170, 190, 220 MeV for a,b and c respectively. Here d,e and f are the Pqgp v/s T 
curves for the same set of baryon chemical potential values. 



In Fig. 8 we show the quark-hadron phase transition diagram. Here we also compare 
our result with the one derived for point like hadrons. It is quite clear that the effect 
of the finite size of baryons is to shift the CEP of the quark-hadron phase transition 
to higher chemical potential values. Also for any baryon-chemical potential value the 
temperature values corresponding to the phase coexistence are seen to drop as compared 
to that for the baryons with no finite size. However one can see the drop is not very 
large, this fact can be directly attributed to the large slope of P-T curves for the 
QGP phase. This makes sure that no matter how large the deviation of P-T curves as 
calculated with RMFTWFS from RMFT is, the intersection point in both cases lies very 
much close together. Now in Table. 1, we list the results of CEP for the quark-hadron 
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100 110 120 130 140 150 160 170 
T (MeV) 

Figure 7: Prmft, Pqgp for a chemical potential /i^ =140 MeV, here given by a and b 
respectively. 'A' is the point of intersection which is just above the P=0 axis. Here 'B' 
represents the point below which the entropy for the hadronic phase is negative and is 
therefore amputated 



phase transition as obtained in different studies and compare them with our findings 



here. Keeping into account the versatility of finite-size scaling 



n 



SS ) [H [13 [IS] and 



that of Re-normalization Group ( RG ) [17] methods in general i\ which allows one to 
have the information about the criticality of a system based solely on a very general 
characteristics and without any detailed description, it is of vital importance to compare 
any model based results with the recently performed FSS analysis of RHIC data [20] . 
In this analysis it was found that the CEP of the quark-hadron phase transition must 
occur only for baryon chemical potential values of fiB > 400 MeV. On comparing with 
our result it is clear that our result for the critical end point C2 lies well above the lower 
bound provided by FSS analysis. This is rather strange in view of the fact that our 
result C2 has been derived in thermodynamic limit of infinite volume and FSS analysis 
of RHIC data has been performed on a system of finite size. This seems to us the direct 



""The former can be derived quite naturally by applying RG techniques to the critical phenomena 
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Figure 8: Quark-hadron phase-diagram, here Ci and C2 are the CEP for hard core 
baryon-radius of rB=Ofm and rB=0.8fm , with co-ordinates Ci= 120, 154 and C2— 
140, 152 respectively. 

consequence of interactions present in the modeL 



Table 1: List of CEP's 



Method 




\/snn 


Tc 


FSS 


~ 400 MeV 


~ 5.75 GeV 




Ci 


380 MeV 




154 MeV 


C2 


420 MeV 




152 MeV 


HRG 


156 MeV 




160 MeV 


QPM-I 


166 MeV 




183 MeV 


QPM-II 


155 MeV 




166 MeV 



One can further clarify the above statement as follows, we know that the size of the 
system comes into picture while calculating the functions of the form, ^ /(p, m) — )■ 
■^V J (PpflpjEj-m), which for the interacting system, under the mean- field approxi- 
mation takes the form, J d^r J d^pf{p,E*,m*). Taylor expanding the function 'f 

9.f, 



we have, V/h^ J d^p[f{p,E,m) — adf/dE — pdf/dm — ...], which can be written as, 
l/h^V{l - ad/dE - Pd/dm - ...) J d^pf{p,E,m) = 1/h^VintJ d^pf{p,E,m). Clearly 
Vint/y < 1- Therefore for a certain specified class of functions f(p,E,m), which have 
well defined n-th derivative and are such that adf/dE, pdf /dm are positive dimension- 
less quantities, a system of interacting particles can be replaced by a non-interacting 
system of reduced volume. 

Now we further compare our results with the analysis done in ref. (21]. Here HRG 
EOS was used for the hadronic phase and bag-model EOS for the QGP phase. It is 
clear that the CEP obtained in this model is much below the lower bound provided by 
the FSS analysis. A similar conclusion can be drawn from the analysis as done in ref. 
[22] where HRG EOS was again used for the hadronic phase and quasi-particle model 
QPM-I and QPM-H [231 El |25l ESI [27] , were used to study the QGP phase. Here again 
the CEP obtained is much below the lower bound provided by the FSS analysis. 
Therefore from the above comparisons it becomes clear that the effect of incorporating 
finite size correction for baryons in the RMFT equation of state has a profound effect 
on the location of critical end point of the quark-hadron phase transition. The addition 
of such effects makes sure that the CEP is well within the bound provided by the FSS 
analysis. More precisely well above the lower bound of FSS analysis. 
Now here it is important to mention that the calculations here have been kept up to 
mean-field level only. This is primarily because of the fact that the couplings are large 
in our model. However this does not negate the fact that the corrections corresponding 
to calculations beyond mean-field level should play its part in the location of the critical 
end point of this quark hadron phase transition. As this requires a detailed study we 
will pick this in future studies. 

VII. Summary 
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In this work we have studied the effect of hard core volumes of hadrons on the quark 
hadron phase transition and in particular on the the critical end point CEP, of such 
a phase transition. The corrections to the hadronic equation of state, for excluded 
volume effect are incorporated in a thermodynamic consistent manner for vander walls 
like interaction. We have found that quantum mechanical generalization of the EOS for 
Boltzmann particles with two particle interactions can be written in a form similar to 
the one obtained for an RMF model with an extra term corresponding to the finite size 
correction. We studied the effect of such hard core volumes of hadrons on the quark 
hadron phase diagram. We obtain that the effect of incorporating the finite hard core 
volumes of the hadrons is to shift the CEP of the quark hadron phase transition to the 
higher chemical potential values. Our findings are consistent with the finite size scaling 
of the RHIC data which negates the CEP with baryon chemical potential values of less 
than 400 MeV. 
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Appendix A.I 



From the symmetry of the pressure 'P' under the transformation, Hb ^ —j^B or ^ 
—/lb, it is clear that the field Uq must have a contribution from a term which is an odd 
function of and n^. An obvious choice would be the function C/Q(ns, n^, T), where 
ub and rib must always appear in product form. However we choose a more general 
form of the function Uq as, 



Uq = C/i(n^,n^^,n^sr) + [/2(n^,nB,n6,n^Sn^«,T) (A.l) 
Uq = Uiin"", T) - U2{n" , ub, Ub, n^'',T) (A.2) 

Now, taking the derivative of pressure P with respect to variable , we get. 



dP _ _ hOUq dPQ 



Also, for the variation of pressure P with variable ub, one has. 



dP dUQ dUQ dPQ 

^ — = -ns^ - rib^ — + -T-^ A.4) 

OUb OUb OHb OUb 



Taking a derivative of A. 3 with hb and A.4 with and by using the independence of 



9.R 



ub, and Ub, □ and subtracting the resulting equations we have, 



Now because of the independence of ra^, n^, rif, the above equation yields, 

d'Uo 



dn^dns 

which after integrating twice yield. 



(A.6) 



Uq = U{n^ ,n^\n^^ ,T) + U{nk,nB,n^\n^^ ,T) + f/(r^•^^ ra^^, T) (A. 7) 

Now using again the independence of pair of variables ^ n^''; , n^^; ra"^^, and 
following the same procedure as above, one obtains. 



Uq = U^{n'',T) + U2{n'\n''^,T) + Us{n,,n^\nB,T) (A.8) 
f/Q = ^7l(n^,T) + f/2(r^^^r^^^,T)-f/3Kr^^^r^B,T) (A.9) 



Appendix B.I 

Let us consider three arbitrary functions A = A{Ax,Ay),B = B{B,j.,By) and C = 
C{Cx,Cy) with the transformation property, Det(|^) 7^ and Det(|^) 7^ 0, then it 
can be shown that for the transformation. A—)- C one has Det(|^) 7^ 0. However if 
^Appendix B.I 
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either Det(||) = or Det(|g) = 0, then it can be shown that, Det(|^) = oM Let 
us consider a transformation A=A(C), and assume that the determinant of Jacobian 



matrices, for this transformation, satisfies. 



(B.i: 



— ^ — y ^ — ^ = (B.2) 

dC^dCy dCydC^ ^ ' 

Now if there exists the transformations, A = A{Bx, By) and B = B{Cx, Cy), then one 
can write, 



dA^ dAy _ dA^ dB^ dAy dB^ ^ dA^^ dB^ dA, dB„ 



dC^dCy dB^dC^dB^dCy dB^dC^dBydCy 



dA^dB^d^dB, dA , . 

dBy dC^ dB^ dCy dBy dC^ dBy dCy ^ ' ' 



dA^ dAy _ dA^ dB^ dAy dB^ ^ dA^ dB,, dAy dBy ^ dA^ dBy dAy dB^ 



dCydC^ dB^dCydB^dC^ dB^dCydBydC^ dBydCydB^dC^ 
d_A^d_B^dAdBy_ 

dBy dCy dBy dC ^ ^ ^ ^ 

therefore, 

dA^dAy dA^dAy fdA.^dAy dA^dAy\ [dB^dBy dB^dBy 



dC^dCy dCydC^ \dB^dBy dBy dBj \dC.^ dCy dCy dCy 

^ :dA\ ^ fdB\ 
--Det { — ] Det i — ] 



dB J \dC J 



-'^"It is for the sake of simplicity that we have kept the functions A,B,C confined to two components 
only. 



Now Eq. B.l becomes, 




(B.5) 



Which imply that either, Det (|^) = or Det (|^) = 0, and also quite simply that, 
Det (1^) = = Det (§§)■ Now if instead of expression B.l, one starts with the 



assumption that for the transformation A=A(C), determinant of Jacobian is non-zero, 
i.e Det 7^ 0, then it is quite clear from B.4 that, Det (||) ^ and Det ^ 0. 
Now since , ub are the independent parameters and riB, n*^^ are independent also, 
therefore the determinant of the jacobian matrices for the transformations relating these 
variables must be non-zero, i.e. 



It therefore follows directly from the above discussion that and n'^^ are also inde- 
pendent, i.e 



Also for the independent set of variables (n^, n^) and (n;,, n^^), it follows from the above 
discussion that {n^,n^^) form another independent pair. 





(B.6) 
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